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18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting # 5: The Fundamental Solution for the Heat Equation 



1. The Fundamental solution 

As we will see, in the case Q = M. n , we will be able to represent general solutions the inhomoge- 
neous heat equation 



n 

(1.0.1) u t — DAu = f, A^^ 2 

i=l 



in terms of /, the initial data, and a single solution that has very special properties. This special 
solution is called the fundamental solution. 

Remark 1.0.1. Note that when Q = R n , there are no finite boundary conditions to worry about. 
However, we do have to worry about "boundary conditions at oo." Roughly speaking, this means 
that we have to assume something about the growth rate of the solution as \x\ — > oo. 



Definition 1.0.1. The fundamental solution Tjj(t, x) to (1.0.1) is defined to be 



(1.0.2) ^(t,x) d = f (4 J t)n/2 e-^, i>0,,Gl», 

where x = f (x 1 , ■ ■ ■ , x n ), \x\ 2 = f Y^i=x( xl ) 2 - 

Let's check that T£>(t,x) solves ( |1 .0. 1 ) when / = in the next lemma. 



Lemma 1.0.1. Tr>(t,x) is a solution to the heat equation (1.0.1) when f = for x G lR n ,t > 0. 



Proof. We compute that d t Y D {t, x) = (- (4i S 3+1 + U7T ^ t)n / 2 jffir) e Also, we compute diV D (t,x) 



-(ld§F2+r e S * and dfT D (t,x) = ( - (4 ^/2+i + 45t (J^W2 + i ) e ^ 
DAT D (t,x) = ( - pjp^ + ^(^Srx)e-^. Lemma 
Here few very important properties of T£)(t,x). 

Lemma 1.0.2. T£>(t,x) has the following properties: 

(1) Ifx^O, then limt^ + ^(t, x) =0 

(2) lim^ Q+ V D (t, 0) = oo 

(3) / RB T D (t, s) d n a; = 1 /or all t > 

Proof. This is a good exercise for you to do on your own. □ 

As we will see, (1) - (3) suggest that at t = 0, Td(0,x) behaves like the "delta distribution 
centered at 0." We'll make sense of this in the next lemma. 

l 



1.0.1 



now easily follows. □ 
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Remark 1.0.2. The delta distribution is sometimes called the "delta function," but it is not a 

function in the usual sense! 

So what is the delta distribution? 

Definition 1.0.2. The delta distribution 6 is an example of a mathematical object called a distri- 
bution. It acts on suitable functions 6(x) as follows: 



1.0.3) 



Remark 1.0.3. The notation (-, ■) is meant to remind you of the L inner product 



1.0.4) 



(f,9) 



f{x)g{x) d n x. 



The next lemma shows that Tu(t,x) behaves like the delta distribution as t — > + . 

Lemma 1.0.3. Suppose that 0(x) is a continuous function on lR n and that there exist constants 
a, b > such that 



(1.0.5) 

Then 



fx) I < ae' 



b\x\ 2 



1.0.6) 



lim 

i-s>0+ 



T D (t,x)(p(x)d n x = 0(0). 



Proof. Using Property (3) of Lemma 1.0.2, we start with the simple inequality 



(1.0.7) 0(0) 



T D (t,x)0(O)d n x 



T D (t,x)(p(x)d n x+ / r D (t,x)(0(O) -<p(x))d n x. 



Let e > be any small positive number, and choose a ball B of radius R centered at such that 
|0(O) — <f)(x)\ < e for x G B (this is possible since is continuous). Then the last term from above 
can be estimated as follows, where B c denotes the complement of B in M. n : 



T D (t,x)(0(O) -<p(x))d n x 
(1.0.8) 

We have thus shown that 



< / r D (t,x)|0(O) -<f)(x)\d n x+ / T D (t,x)\(f)(0) - (j)(x)\d n x 

B JB C 



< / T D (t,x)ed n x + 



B 



T D (t,x)d n x+ / T D (t,x)\4>(x)\d n x 



B' 



< e+ |0(O)| / T D (t,x)d n x+ [ r D (t,x)\<f)(x)\d n x. 

JB C Jb c 



(1.0.9) 0(0) 



T D (t,x)(t)(x)d n x 



< e + 



T D (t,x)d n x+ / Y D {t } x)\(t)(x)\d n x. 

B c JB C 



To estimate the final term on the right-hand side of (1.0.8), we take advantage of the spherical 

dof 



symmetry of T(t, x) in x. More precisely, we introduce the radial variable r = \x\ = \ZY^= 



x 
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and recall from vector calculus that for spherically symmetric functions, d n x = C n r n ~ l dr where 
C n > is a constant. Therefore, using the assumed bound \<j)(x)\ < ae br , we have that 



1.0.10) 

T D (t, x)\(j){x) | d n x < C' n t~ n/2 I r n ~ l e-^- 



r=R 



-n/2 



p n - l e~ p dp, 



where C' n > and C" > are constants. To deduce the second inequality in (1.0.10), we have 



made the change of variables r = p( 



J_ _ 6 )-l/2 



4Dt 



P t 



i/2rJ_ _ 



bt) 1/2 . Now since RJ jj^ - b ->■ oo 



as t — > + , it easily follows from the last expression in (1.0.10) that j BC T£,(t,x) d n x goes to as 
/ >()'. 



The second term on the right-hand side of (1.0.8) can similarly be shown to go to as t — > + . 



Combining the above arguments, we have thus shown that for any e > 0, 



(1.0.11) 

We therefore conclude that 



lim sup 



0(0)- / T D (t,x)(p(x)d n x 



< e. 



(1.0.12) 
as desired. 



lim 

t->o+ 



0(0) 



r D (t, x)4>(x) d n x 



□ 



Remark 1.0.4. Lemma fl.0.3l can be restated as 



1.0.13) 



lim(r D (t,-),0(-)) = (5(-),0(-)) = 0(O). 



On the left, (, ) means the integral inner product, whereas in the middle it has the meaning of 
(1.0.3). We sometimes restate (1.0.13) as 



(1.0.14) 

Let's summarize the above results. 



lim T D (t, x) = 6(x). 



Proposition 1.0.4 (Properties of Tj)(t, x)). Tjj(t,x) is a solution to the heat equation (1.0.1) 
( with f = 0) verifying the initial conditions 



1.0.15) 



lim T D (t, x) = 8(x). 



1.1. Solving the global Cauchy problem when n — 1. Let's see how we can use Td to solve 
the following initial value (aka Cauchy) problem: 



(1.1.1) u t -Du xx = 0, (t, x) G (0, oo) x K, 

u(0,x) = g(x). 

We will make use of an important mathematical operation called convolution. 



4 
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Definition 1.1.1. If / and g are two functions on IR n , then we define their convolution / * g to be 
the following function on M. n : 



(1.1.2) 



(f*g)(x) 



dcf 



f(y)g(x -y)d n y. 



Convolution is an averaging process, in which the function f(x) is replaced by the "average value" 
of f(x) relative to the "profile" function g(x). 

The convolution operator plays a very important role in many areas of mathematics. Here are 



two key properties. First, by making the change of variables z = x — y, d n z = d n y in (1.1.2), we see 
that 



1.1.3) 



f{y)g{x - y) d n y 



f(x-z)g(z)d n z=(g*f)(x), 



which implies that convolution is a commutative operation. Next, Fubini's theorem can be used to 
show that 



(1.1.4) 

so that * is also associative. 



* (9 * h) = (/ * g) * h, 



Remark 1.1.1. According to (1.0.3) and (1.1.3) 



(1-1.5) (/ * 6)(x) = (6{y), f(x - y)) y = f(x), 

so that in the context of convolutions, the 5 distribution plays the role of an "identity element." 

The next proposition is a standard fact from analysis. It allows us to differentiate under integrals 
under certain assumptions. We will use it in the proof of the next theorem. 

Proposition 1.1.1 (Differentiating under the integral). Let I(a,b) be a function on 1 x R. 
Assume that 



(1.1.6) / \I(a,b)\da < oo 

Jr 

for all b belonging to a neighborhood of b and define 



(1.1.7) JWSflMCa. 

Jr 

Assume that there exists a neighborhood Af of b such that for almost ev en£] a, dbl(a, b) exists for 
b G N '. In addition, assume that there exists as function U(a) (defined for almost all a) such that 
for b G Af, we have that \dbl(a, b)\ < U(a) and such that 

In a measure theory course, you would learn a precise technical definition of "almost every." For the purposes of 
this course, it suffices to know the following fact: if a statement holds for all a except for those values of a belonging 
to a countable set, then the statement holds for almost every a. The main point is that the function I(a, b) does not 
have to be "well-behaved" at every single value of a; it can have some "bad a spots," just not too many of them. 
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(1.1.8) / U(a)da < oo. 

Jr 

Then near J{b) is differentiable near b , and 



d b J{b) = / d b I(a, b) da. 



(1.1.9) 



Remark 1.1.2. An analogous proposition is true for functions I(a,b) defined on 



x 



Theorem 1.1 (Solving the global Cauchy problem via the fundamental solution). Assume 
that g(x) is a continuous function on W 1 that verifies the bounds \g(x)\ < ae b ^ , where a, b > are 
constants. Then there exists a solution u(t, x) to the homogeneous heat equation 

(1.1.10) u t -DAu = 0, (t>0,xeK n ), 

u(0,x)=g(x), xeW 1 

existing for (t,x) G [0, T) x IR n , where 



(1.1.11) 

Furthermore, u(t,x) can be represented as 



T 



def 1 



ADb 



1.1.12) 



u{t,x) = \g(-)*r D (t,-)](x) 



9(y)^D(t,x- y) d n y 



(AnDt) n l 2 



g(y)e d n y. 



The solution u(t,x) is of regularity C°°((0, jj^A x]R n ) (i.e., it is infinitely differentiable). Finally, 
for each compact subinterval [0,T'] C [0, T), there exist constants A,B>0 (depending on the 
compact subinterval) such that 



1.1.13) 



\u{t,x)\ < Ae B ^ 2 



for all (t,x) G [0, T'] x ]R n . The solution u(t,x) is the unique solution in the class of functions 



verifying a bound of the form (1.1.13) 



Remark 1.1.3. Note the very important smoothing property of diffusion: the solution to the 
heat equation on all of M. n is smooth even if the data are merely continuous. 



Remark 1.1.4. The formula (1.1.12) shows that solutions to (1.1.10) propagate with infinite 



speed: even if the initial data g(x) have support that is contained within some compact region, 



(1.1.12) shows that at any time t > 0, the solution u(t,x) has "spread out over the entire space 



MV In contrast, as we will see later in the course, some important PDEs have finite speeds of 
propagation (for example, the wave equation). 



6 
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Proof. For simplicity, we only give the proof in the case n = 1. The basic strategy of the proof is to 
analyze the behavior of Tn{t,y) in detail. 

Let u(t, x) be the function defined by ( jl.l. 12 ) . The argument that follows will show that the right- 
hand side of ( jl.l. 12 ) is finite (and more). In fact, let us first demonstrate the bound (1.1.13). To this 
end, let e > be any positive number. Then using the simple algebraic estimate \2xy\ < e~ 1 x 2 + ey 2 , 
we deduce the inequality 



(1.1.14) 



\x - y\ 



x 



2xy + y 2 < (1 + e^)x 2 + (1 + e)y 2 . 



Using (1.1.14) and the assumed bound on |<?(-)|, we deduce that 



(1.1.15) 



Ax -y)\< ae b \ x ~ y \ 2 < ae^ 1 ^ V 1 ^^ 



Using (1.1.15) and the fact that L n g(y)Tc{t,x — y) dy = J Rn g(x — y)T r>{t,y) dy (i.e., that convo- 
lution is commutative), we have the following estimates: 



1.1.16) 



u(t,x)\< / Ux-y^Tnit^dyKae^- 1 ^ 2 / e^ b ^T D (t, y) dy 



dy 



ae 



(l+e- 1 )^! 2 



f 1 r 1/2 e~[^~ (1+e)6 ] : ' 



<oo 



where A > is an e— dependent constant, and in the next-to-last step, we have made the change of 
variables z = [^jyj. — (1 + e)b] l ^ 2 y = t~ x l 2 — + e )^] 1 ^V- Note that this change of variables 

Since e is allowed to be arbitrarily small, we have thus 



i 



is valid as long as < t < 
demonstrated an estimate of the form (|1. 1.131). 



47rD(l+e)6 



Let's now check that the function u(t,x) defined by (1.1.12) is a solution to the heat equation 

and also that it takes on the initial conditions g(x). To this end, let £ = f dt — Dd 2 . We want to show 
that Cu(t,x) = for t > 0, x G lR n and that u(t, x) — > g(x) as t \, 0. Recall that by Proposition 



1.0.4, CT D (t, x) = for t > 0, x E R. For t > 0, x G R, we have that 



1.1.17) 



Cu(t, x) 







g(y) £T D (t,x- y) dy = 0. 



To derive (1.1.17), we have used Proposition 1.1.1 to differentiate under the integral; because of 



rapid exponential decay of Td(-, ■) in its second argument as the argument goes to oo, one can use 
arguments similar to those given in the beginning of this proof to check that the hypotheses of the 
proposition are verified. 

Similarly, the fact that u G C°° ((0, x R ) can be derived by repeatedly differentiating with 



respect to t and x under the integral in (1.1.12). 
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Furthermore by (1.0.15) and (1.1.5), we have that 
(1.1.18) 



lim u(t,x) = lim (#(•) * T D (t, -))(x) = (g * 6)(x) = g(x). 
t->o+ t^o+ 



The question of uniqueness in the class of solutions verifying a bound of the form (1.1.13) is 



challenging and will not be addressed here. Instead, with the help of the weak maximum principle, 
you will prove a weakened version of the uniqueness result in your homework. 

□ 



In the next theorem, we extend the results of Theorem LI to allow for an inhomogeneous term 
f(t,x). 



1.1 



Also assume 



Theorem 1.2 (Duhamel's principle). Let g(x) andT = be as in Theorem 
that f(t, x), dif(t, x), and d{djf(t, x) are continuous, bounded functions on [0, T) x W L for 1 < i, j < 
n. Then there exists a unique solution u(t,x) to the inhomogeneous heat equation 



(1.1.19) ut- DAu = f{t,x), (t,x) G (0, oo) x R, 

u(0, x) = g(x), i6l 

existing for (t,x) G [0, T) x R. Furthermore, u(t,x) can be represented as 



(1.1.20) u(t,x) = (r D (t,.)*g)(x)+ / (r D (t-s,-)*f(s,-))(x)ds. 

Jo 

The solution has the following regularity properties: u G C°([0,T) xl)(l C 1,2 ((0,T) x R). 

Proof. A slightly less technical version of this theorem is one of your homework exercises. □ 

2. Deriving T D (t, x) 

Let's backtrack a bit and discuss how one could derive the fundamental solution to the heat 
equation 



(2.0.21) d t u(t,x) — DA x u(t,x) — 0, (t, x) G [0, oo) x E n . 

As we will see, the fundamental solution is connected to some important invariance properties 



associated to solutions of (2.0.21). These properties are addressed in the next lemma. 



Lemma 2.0.2 (Invariance of solutions to the heat equation under translations and par- 



abolic dilations). Suppose that u(t,x) is a solution to the heat equation (2.0.21). Let A, to G M be 

constants, and xq G M n . Then the amplified and translated function 



(2.0.22) 



u* (t, x) = Au(t — t , x — Xo) 



is also a solution to (2.0.21). 



Similarly, if A > is a constant, then the amplified, parabolically scaled function 



(2.0.23) 

is also a solution. 



u*(t,x) = f Au(X 2 t,Xx) 



cS 
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Proof. We address only the case (2.0.23), and leave (2.0.22) as a simple exercise. Using the chain 



rule, we calculate that if u is a solution to (2.0.21), then 



(2.0.24) d t u*(t,x)-Au*(t,x) = X 2 A{{d t u)(X 2 t,Xx) - (DAu)(X 2 t, Ax)} = 0. 

Thus, u* is also a solution. 



□ 



We would now like to choose the constant A in (2.0.23) so that the "total thermal energy" of u* 
is equal to the "total thermal energy of" of u. 

Definition 2.0.2. We define the total thermal energy T{t) at time t associated to u(t, x) by 



(2.0.25) 



T{t) 



def 



u(t, x) d n x. 



It is important to note that for rapidly-spatially decaying solutions to the heat equation, T{t) is 
constant. 

Lemma 2.0.3. Let u(t,x) G C 1,2 ([0,oo) x IR n ) be a solution to the heat equation —d t u(t,x) + 
Au(t, x) = 0. Assume that at each fixed t, lim^i^oo |x| n_1 |V x M(t, x)\ = 0, uniformly in x. Further- 
more, assume that there exists a function f(x) > 0, not depending on t, such that \d t u\ < f(x) and 
such that L n f(x) d n x < oo. Then the total thermal energy of u(t,x) is constant in time: 



(2.0.26) 



T(t) = T(0). 



Proof. Let T(t) = f j Rn u(t, x) d n x denote the total thermal energy at time t. The hypotheses on 
ensure that we can differentiate under the integral and use the heat equation: 



(2.0.27) 



d_ 
dt 



T(t) 



d t u(t, x) d n x 



Au(t, x) d n x = lim 



Au(t, x) d n x, 



Sfl(O) 



where -Br(O) C M n denotes the ball of radius R centered at the origin. Then with the help of 
the divergence theorem, and recalling that da = R^^du along <9I?r(0), where u denotes angular 
coordinates along the unit sphere dBi(0), we conclude that 



(2.0.28) lim / Au(t,x)d n x= lim / Vj^u(t,a)da 

R ^°° JB R (0) R ^°° JdB R (0) 

= lim / R n - 1 Vfiu{t,Ru)du 



I 



lim R W fju{t } Rjj) doj 



aBi(o) 



R— >oo 



0du = 0. 



aBi(o) 



In the last steps, we have used the following basic fact from analysis: the condition limia-i. 
uniformly in u allows us to interchange the order of the limit and the integral. 



\x 



n— 1 



V x ii(t, x) 
□ 



We now return to the issue of choosing constant A in (2.0.23) so that the total thermal energy of 



u* is equal to the total thermal energy of of u. Using the change of variables z = Ax, and recalling 
from multi-variable calculus that d n z = ~h n d n x 1 we compute that 
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(2.0.29) / u*(t,x)d n x = A u(D 2 X 2 t,Xx) d n x = AX' 12 u(X 2 t,z)d n z. 

Observe that that L„ u(X 2 t, z) d n z is in fact the mass of u. Thus, we choose A = X n , which results 



in 



(2.0.30) 



u*(t,x) = X n u(D X t, Xx). 



Motivated by the parabolic scaling result (2.0.23|), we now introduce the dimensionless variable 
(2.0.31) C 



def X 



where we have used the fact that the constant D has the dimensions of [length 2 ] /[time]. Note that 
( is invariant under the parabolic scaling t — > X 2 t, x — > Xx. 

We now proceed to derive the fundamental solution. For simplicity, we only consider the case of 
1 + 1 spacetime dimensions. We will look for a fundamental solution of the form 



(2.0.32) 



r D (t,x) 



v(0, 



where V(() is a function that we hope to determine. Admittedly, it is not easy to completely 
motivate the fact that r^(t,x) should look like (2.0.32). We first note that since we would like to 
achieve f R T D (t,x) = 1, the change of variables (2.0.31) leads to the following identity: 



(2.0.33) 



r D (t,x) 



' v(-?L)<l.r 



v(C) d(. 



Dt WDt' 

Next, since T D (t,x) is assumed to solve the heat equation, we calculate that 

(2.0.34) = d t T - AT = - j nc) + + lv(C)}. 

Therefore, V must be a solution to the following ODE: 



no + \cv(o + \v(o = o. 



(2.0.35) 

Since we want ro(t,x) to behave like the 5 distribution (at least for small t > 0), we demand that 



(2.0.36) V{() > 0. 

Furthermore, since we want Tu(t, x) to rapidly decay as |x| — > oo, we demand that 



(2.0.37) 



V(±oo) 



0. 



We also expect that ideally, ^(C) should be an even function. Furthermore, it is easy to see that if 



V(() is a solution to (2.0.35), then so is W(() = f V{— (). Thus, it is reasonable to look for an even 
solution. Now for any differentiable even function V((), it necessarily follows that V'(0) = 0. Thus, 
we demand that 



10 
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(2.0.38) 



7(0) = 0. 



We now note that (2.0.35) can be written in the form 

d 



(2.0.39) 



which implies that V'(Q + |£7(C) is constant. By setting ( = in and using (2.0.38), we see that 
this constant is : 



(2.0.40) 



HO + g WO = o. 



Now the first-order ODE (2.0.40) can be written in the form 

d 



(2.0.41) 

which can be easily integrated as follows: 



hv(0 — i<, 



(2.0.42) 



In 



\V(0) J 4 

(2.0.43) 7(C) = 7(0)e-K 2 . 

To find 7(0), we use the relation (2.0.33), and the integral identity^] 



(2.0.44) 1= I V(0)e-^ 2 d( = 27(0) / e~ a2 da = 27(0)v^. 

Jr <= 2a Jr 

Therefore, 7(0) = and 

1 



(2.0.45) 



Arc 



e 4- 



Finally, from ( |2.0.32[ ) and (J2XL45J), we deduce that 

T D (t,x) = - 



(2.0.46) 
as desired 



1 _^ 

e 4t 



Ant 



2 Let I = Le x * dx. Then I 2 = j' f R e. ^ dx dy , and by switching to polar coordinates, we have that 
I 2 = 2tt J™ q re- r2 dr = tt. Thus, I = ■ s pn. 
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